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Abstract: Phantom and mouse experiments of time-domain fluorescence 
tomography were conducted to demonstrate the total light approach which 
was previously proposed by authors. The total light approach reduces the 
computation time to solve the forward model for light propagation. Time- 
resolved temporal profiles were acquired for cylindrical phantoms having 
single or double targets containing indocyanine green (ICG) solutions. The 
reconstructed images of ICG concentration reflected the true distributions of 
ICG concentration with a spatial resolution of about 10 mm. In vivo 
experiments were conducted using a mouse in which an ICG capsule was 
embedded beneath the skin in the abdomen. The reconstructed image of the 
ICG concentration again reflected the true distribution of ICG although 
artifacts due to autofluorescence appeared in the vicinity of the skin. The 
effectiveness of the total light approach was demonstrated by the phantom 
and mouse experiments. 
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1. Introduction 

Over the last decade, a large number of research and development has conducted for the 
technology of molecular imaging, and one of the modalities of molecular imaging is 
fluorescence diffuse optical tomography (FDOT) [1,2]. FDOT is an emerging technology to 
image spatial distributions of fluorescence parameters in tissues which are irradiated by a 
near-infrared (NIR) excitation light source and emit fluorescence light [3,4]. The fluorescence 
parameters include fluorophore concentration, fluorescence quantum yield and lifetime. 
Accumulation into and secretion from specified organs of the target drug labeled with a 
fluorophore can be investigated by tracing the temporal change of the fluorophore 
concentration. Fluorescence quantum yield and lifetime can give information about the tissue 
microenvironments such as temperature, viscosity, pH, etc. on which they are dependent. 

Essentially, FDOT is an extension of the technology of diffuse optical tomography (DOT) 
[5,6]. DOT reconstructs the tomographic images of absorption and scattering characteristics 
inside the tissue from the data of detected light intensities which are irradiated onto and 
reemitted from the tissue surface using the reconstruction algorithm based on the inversion 
technique. Usually, near-infrared light is used due to its high transmission capability through 
biological tissues, and the images of absorption characteristics are converted to the images of 
physiological information through the absorption spectra of chromophores in biological 
tissues, such as hemoglobin and myoglobin. The technique of DOT can be classified into 
three modes according to the type of source light, (i) continuous wave (CW), (ii) frequency 
domain and (iii) time domain. Because the number of useful data for reconstruction increases 
in this order, time domain mode is believed to give the images with highest quality with the 
costs of sophisticated instruments for measurement and long computation times for 
reconstruction [7]. 

These three modes can be used for FDOT as well. So far, CW and frequency domain 
modes have been mainly used for FDOT and many results have been reported [8-17] because 
time domain mode needs complicated and expensive instruments and takes longer time to 
acquire data for reconstruction. However, it is difficult for CW mode [8-10] to measure the 
fluorescence lifetime which has the order of nanosecond. Frequency and time domain modes 
can reconstruct both the fluorescence concentration and lifetime, but so far many FDOT 
investigations have been conducted using frequency domain mode [11-17] while those using 
time domain mode are limited. FDOT investigations using time domain mode have been 
conducted by simulations [18-21] first, and some experimental works have been reported 
recently [22-25]. 

In reconstruction for FDOT, the forward model of light propagation is based on the photon 
diffusion equations for excitation and emission light. Particularly for time domain mode, the 
emission source term in the photon diffusion equation for emission light is a convolution 
integral of the excitation light fluence and the fluorescence lifetime function. This integral 
makes the computation of the forward model complex and long. To avoid the computation of 
the convolution integral, the total light approach was proposed previously and simulation 
results were presented [21]. 

The purpose of this study is to experimentally demonstrate the total light approach for 
time domain FDOT [21] by showing the experimental results using phantoms and mouse. 
Time-resolved measurement employing a time correlated single photon counting system 
requires longer time for data acquisition with a satisfactory S/N ratio than CW or frequency 
domain modes. But it has been reported that high quality images of DOT can be reconstructed 
with the use of the time-resolved data at the costs of long computation and large memory [7], 
and time domain FDOT can provide the images of the fluorescence characteristics such as the 
fluorophore concentrations, the fluorescence quantum yield and lifetime [20,25]. In the 
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previously proposed total light approach, however, the fluorescence quantum yield and 
lifetime are assumed to be known and only the fluorophore concentrations were reconstructed. 
Therefore, the total light approach does not fully utilize the advantages of time domain mode. 
But the previous report of FDOT simulations employing the total light approach easily 
reconstructed the fluorophore concentration with the measured mean time-of-flight data as the 
featured input data for reconstruction. To our knowledge, none of the reported experimental 
investigations using time-resolved measurement [22-25] has used the mean time-of-flight 
data obtained from the measured time-resolved profiles. Niedre et al. [22] used the early 
photons in the time period of 100 ps for reconstruction, Northdurft et al. [23] and Kepshire et 
al. [24] converted the measured time domain data to frequency domain data for 
reconstruction, and Gao et al. [25] employed the generalized pulse spectrum technique for 
reconstruction where the measured time -resolved data were Laplace transformed. As the first 
step toward the full use of the time -resolved measurement data for FDOT, this study reports 
the experimental results of time domain FDOT to demonstrate the total light approach using 
the measured mean time-of-flight data as the featured input data. 

In the demonstration experiments in this study, we focus our attention only on imaging the 
fluorophore concentrations within the objects of cylindrical phantoms and a mouse in vivo. 
Indocyanine green (ICG) was used as a fluorophore, and capsules containing ICG were 
inserted inside the phantoms and mouse abdomen. The data were acquired with a time- 
resolved measurement technique, and the tomographic images of ICG concentrations in the 
capsules were reconstructed by use of the total light approach. In the following, the total light 
approach is briefly explained, the experimental method including the instruments and objects 
are described, the experimental results are reported, and finally the conclusions are stated. 

2. Image reconstruction using the total light approach 

2. 1 Formulation of the forward problem by total light 

Image reconstruction is performed by solving the inverse problem where the forward problem 
for propagation of fluorescence excitation and emission light must be modeled. Light 
propagation in random media is strictly described by the radiative transfer equation (RTE), 
and some investigations of FDOT based on RTE have been reported [26,27] but limited to 
CW mode probably due to the heavy computational load. Then, the photon diffusion equation 
which is the diffusion approximation of RTE is frequently used for the model of light 
propagation in biological objects with diameters larger than 10 mm [5,6]. The time dependent 
photon diffusion equations for excitation and emission light are expressed by Eqs. (1) and (2) 
[11,28]. 



-V • (Z)(r) V) + ju a (r) + eN(r) + 



1 d_ 
c dt 



(1) 



-V-(D(r)V) + // fl (r) + 



1_ d_ 
c dt 



0 m (r,t) = ] €N(r)\y x (r,t')-- e xp{-(t-t')/T}dt', (2) 



where r [mm] is the position vector, t [ps] the time, 4>(r,?) [W/mm ] the fluence rate, q{v,t) 
[W/mm 3 ] the internal light source, c [mm/ps] the speed of light in the object, Z)(r) [mm] the 
photon diffusion coefficient, /i a (r) [mm -1 ] the absorption coefficient, e [mnf'uM" 1 ] the 
specific extinction coefficient of the fluorophore for the excitation light, N(r) [uM] the 
concentration of the fluorophore, r [ns] the lifetime of the fluorophore, y the quantum yield of 
the fluorophore, subscripts x and m denote the excitation and emission light. The photon 
diffusion coefficient is expressed as Z)(r) = l/(3,«/(r)) using the reduced scattering coefficient, 
fi s '(r) [29]. Because the wavelength of excitation is different from that of emission Z)(r) and 
fi a { r ) m Eq. (1) af e different from those in Eq. (2). But the difference in the wavelength is as 
small as about 10 nm, Z)(r) and fi a {r) for excitation wavelength are assumed equal to those for 
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emission wavelength. Also, the specific extinction coefficient, e, and the quantum yield, y, are 
assumed to be known and homogeneous within the object. 

The right hand side of Eq. (2) is the emission term which is the convolution of the fluence 
rate of excitation light, <£> x (r,t), and the lifetime function, exp(-tk)/z, and needs long 
computation time. Now the total light approach is introduced as follows [21]. Assuming an 
ideal fluorescence with infinitesimal lifetime, i.e. z = 0, modifies Eq. (2) as Eq. (3), 



-y.(D(r)V) + Ma (r) + ~ 
c ot 



Ol(r,0 = ^(r)O x (r,0, (3) 



where <f>* m (r,t) is the fluence rate of emission when x = 0, namely the zero-lifetime emission, 
and is related to the real fluence rate of emission, O m (iy), as Eq. (4), 



O m (r, 0 = f o;„ (r, 0 • - exp {-(* - 1') I r}dt'. (4) 



T 

Now we define the total light O^r,?) as the following Eq. (5), 

O r (r,0 = O,(r,0 + -Ol(r,0- (5) 

r 

Then the photon diffusion equation for the total light, O^r,?), is derived from Eqs. (1) and (3) 
as the following, 



c at 



0 T (r,t) = q x (r,t). (6) 



Physically, the total light is understood as the sum of excitation and emission light under the 
extreme fluorescence condition where all the absorbed energy of excitation light is converted 
to fluorescence light (y = 1) with the zero-lifetime (r = 0). Introducing the total light converts 
the combined fundamental equations from Eqs. (1) and (2) to Eqs. (1) and (6). The difference 
between Eqs. (1) and (6) is only the term of absorption of excitation light, sN(r), and thus the 
solutions of <t> x (r,?) and <£>j(r,t) from Eqs. (1) and (6) are numerically obtained more easily and 
faster than the solutions of <& x (r,t) and 4> m (r,?) from Eqs. (1) and (2) using the finite element 
method (FEM) [30]. The fluorophore concentration, N(r), is obtained by reconstructing fi a {r) 
+ e N(r) in Eq. (1) and (i a ( r ) m Eq. (6) first, then by dividing the subtraction {[ fi a (r) + e A^r)] 
-fi a (r)} bye. 

The boundary conditions for Eqs. (1), (2) and (6) are given as the following, 

-D(r 4 )V,® v (r t ,0 = ^* v (r 4 ,0, (y = x,m,T), (7) 

where is the position at the boundary, V„ the gradient outward normal to the boundary, 
A = (\ + R f )l (\-R f ) the factor with the internal reflectivity at the boundary, Rf, and the 
subscript v represents x, m, or T. R f is given by Eq. (8), 

R f =-1.440«~ 2 +0.710« 1 + 0.668 + 0.0636«, (8) 

where n is the ratio of the refractive index of the object to that of the environment. 
The initial condition is given by Eq. (9), 

<D„(r,0 = 0 for?<0, (y = x,m,T). (9) 

O x (r,?) and O^r,?) are calculated from Eqs. (1) and (6) under the boundary and initial 
conditions (7) and (9). Then ®,„(r,0 is calculated from Eqs. (5) and (4) as the solution of the 
forward problem. 
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The measured powers of the excitation, fluorescence and total light, TJr b ,i), T„,(r b ,t), 
Tj{r h ,t), are the fluxes of the fluence rates of the excitation, fluorescence and total light, 
<£> x (r,t), <E>,„(iy) and O^r,?), at the boundary as expressed by Eq. (10), 

T v (r h ,t) = -D(r h )V^ r (r h ,t), (v = x,m,T). (10) 

From Eq. (7), the measured powers, r v (rt»0s become proportional to the fluence rates, 
4>,,(iy), as shown in Eq. (1 1), 



r„(r 4 ,o =— ®,(r 4 ,0. (y = x,m,T). 



(11) 



2. 2 Image reconstruction 



The tomographic images of the fluorophore concentrations are obtained by solving the inverse 
problems [5] reconstructing fi a {r) + e N(r) in Eq. (1) and nJj) in Eq. (6) from the input data of 
the featured data extracted from the temporal profiles of the excitation and fluorescence light 
measured by a time-resolved measurement system. Figure 1 shows the concept of the process 
of image reconstruction in this study. Firstly, the powers of the excitation and fluorescence 
lights measured at the multiple positions at the object surface, r v (r b ,?) and T,„(r b ,f), are 
obtained by time-resolved measurements. Secondly, a deconvolution process using the 
Wiener filtering [31] calculates the power of the zero-lifetime emission, r*,„(rb,0, by use of 
Eqs. (4) and (11) as well as the known lifetime, r. This deconvolution and noise suppression 
process is carried out for only the measured time-resolved data and the process is done prior 
to iteration for reconstruction taking only a few second in our image reconstruction process. If 
the total light approach is not used the convolution of Eq. (4) for the forward calculation 
should be done for every FEM node and every iterative updating process. Therefore, the total 
light approach reduces the calculation load for reconstruction. Then the measured total light, 
Tj{r b ,t), is obtained using Eqs. (5) and (1 1). It is known that full use of the temporal profiles in 
the reconstruction process improves the quality of the reconstructed images at the cost of long 
computation time. Here we employ the mean times of flight of the excitation and total light 
defined by Eq. (12), <t(r b )> x and <t(r b )> T , as the input data for image reconstruction as the 
previous paper which proposed the total light approach [21], 



rtT v (r b ,t)dt 

<f(»*)>v=7= . {v=x,n 

J 0 KMdt 



(12) 



Measurement 



Inverse calculation 




Fig. 1 . Concept of fluorescence tomography. It consists of measurements of the excitation and 
fluorescence emission light and image reconstruction process by solving the inversion problem. 
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The inverse calculation is performed as the following. Firstly, the photon diffusion 
equations of light propagation are numerically solved with assuming the distributions of the 
optical properties, fi ax {r) ( = ju a (r) + eN(r)) and /n a i(r) ( = fi a ( r )), using FEM. From the 

calculated powers of the excitation and total light, f x (r h ,t) and f T (r b ,t), the calculated mean 
times-of- flight, <t(r b )> x and <t(r b )> T are obtained. Although the diffusion coefficient, 

D(r), varies with organs it is given homogeneously as 0.39 mm as a known value in this study 
by assuming that the differences in the scattering coefficients of biological tissues are not 
large. This assumption can be one of the causes of errors in the reconstructed images. Then 
the measured mean times-of-flight of the excitation and total light, <t{r h )> x and <t{r h )> T , are 

compared with the calculated ones, <i(r h )> x and <t(r b )> T , and the distributions of the 

optical properties, fi ax (r) and fi aI (r), are updated using Eq. (13), 



where M v and F v (%0 are vectors with their elements consisting of < t(r b ) > v and < i(r h ) > v , 
respectively, Xk is a vector with its elements consisting of fi ax (r) and /u a i(r) of the finite 
element in FEM at the k-th iteration, S/t is the update of Xh J v (5Ct) is the Jacobian matrix 

which describes the changes in <t(r b )> v when fi av (r) of each finite element changes 
infinitesimally. The update Sxk is calculated with the algebraic reconstruction technique 



where the superscript, j, denotes the y'-th measurement position with the maximum of the 
product of the source and detector numbers, 3' v J) (% k ) is the y'-th row vector of the matrix 
J v {ik ) > ^ W is the y'-th element of the vector M„ - F„ {% k ) . 

The solutions of fi ax (r) and fi a i(r) are obtained when the iteration reaches convergence 
where the updates of /u ax (r) and /u a i(r) becomes less than a specified small number. From the 
solutions of fiax(r) = fi a (r) + e N(r) and jU aI (r) = fi a (r), the tomographic image of N(r) = (ju ax (r) 
- fi a (r))f s is finally obtained. 

The advantage of the total light approach is to reduce the computational time for Eq. (4) as 
mentioned above. It takes a long computation time to calculate the convolution in each 
updating process in many iterations, while only one calculation of deconvolution prior to the 
iterative process is required by the total light approach. In this study, we use only the time-of- 
flight data for reconstruction. But the total light approach allows the users to use the entire 
time-resolved data with relatively small computation load. The reconstruction procedure 
presented in this paper is basically the double DOTs for the excitation and total lights. Any 
kind of the featured data used in several DOT schemes can be applied to reconstruct the 
fluorophore concentrations. Time-resolved profile is affected by the absorption coefficient. So 
the features of the time-resolved data reflect the absolute value of the optical properties. 
Therefore DOT can reconstruct the absorption coefficient by using the featured data. Mean 
time-of-flight data is not dependent on the intensity of the detected light. Hence the 
calibrations of source and detector are not so important. A possible disadvantage of the 
proposed method is that the method inherits low spatial resolution of DOT. 

Reconstruction of the fluorophore concentrations only can be done with CW systems. 
Time-resolved systems are not necessarily needed for reconstruction only the fuorophore 
concentration. Therefore, the time-resolved measurement system employed in this study is 
over-specified for reconstruction of only the fluorophore concentration. The fluorescence 
lifetime can be reconstructed by better use of the measured time-resolved data. Simultaneous 
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M v -F v ( lk ) = J v ( lk )S lk (v = x,T), 



(13) 



(ART) as Eq. (14), 




(14) 



reconstruction of the fluorescence yield (i.e., fluorophore concentration) and lifetime is 
possible [25], but the reconstruction of the lifetime will be more easily performed if the 
distribution of the fluorophore concentration is known as prior information. Now, the 
fluorophore concentration reconstructed first by the total light approach can provide prior 
information for reconstruction of the lifetime using Eq. (2). This process may be another 
advantage of the total light approach applied to time-domain FDOT. 

3. Experimental setup and objects 

3.1 Fluorescence characteristics of fluorophore 

Indocyanine green (ICG, Sigma- Aldrich) was employed as the fluorophore in this study. ICG 
is approved for injection to humans and widely used clinically for investigating hepatic 
function and retinal angiography, etc. It has absorption and emission spectra in the near- 
infrared wavelength range where light is more transparent to biological tissues. This property 
is very important for fluorescence tomography to be applied to in vivo animals. Figure 2 
shows the measured absorption and emission spectra of ICG contained in 0.8% Intralipid 
solution which has scattering characteristics similar to biological tissues. The peak 
wavelengths of absorption and emission are 767 nm and 820 nm, respectively. The specific 
extinction coefficient, fluorescence lifetime and quantum yield around the wavelength of 800 
nm are given as e = 0.04 mm _1 /uM, r = 0.47 ns, y = 0.03 from the literatures [32-34]. 
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Fig. 2. Absorption and emission spectra of ICG mixed in Intralipid. 



3.2 Time-resolved measurement system 

Figure 3 shows the schematics of the picosecond time-resolved measurement system used in 
this study. It had a light source of a diode laser emitting ultrashort pulse light with a 
wavelength of 759 nm, a pulse width of about 100 ps, a repetition rate of 5 MHz and an 
average power of 0.25 mW. The temporal profiles were measured by a time-correlated single 
photon counting (TCSPC) system which had a time-resolution of about 100 ps. Please refer to 
a literature for the details [35]. Multiple source and detector fiber bundles were used for 
phantom and mouse experiments, respectively. 12 source and detector bundles were 
alternatively attached onto the surface of the cylindrical phantom in the plane perpendicular to 
the cylinder axis with an equal spacing. For in vivo mouse experiments, 8 source and detector 
fiber bundles were attached to the surface of a mouse holder as well. The outer diameters of 
both the source and detector fiber bundles were 4.0 mm. 

The ultrashort pulse light from the light source was introduced into one source fiber 
bundle, and the lights collected by the detector fiber bundles were led to the TCSPC system 
one by one through a switching device. By switching the source fiber bundle one by one, 12 * 
12 (phantom experiments) or 8 x 8 (mouse experiments) temporal profiles were acquired for 
each of excitation and emission lights. When measuring emission light a long pass filter 
(IR82, Fuji film) was used to cut excitation light and a total of 12 x 12*2 or 8*8x2 
temporal profiles were obtained for one image. Figure 4 shows the transmission spectrum of 
the filter as well as the excitation wavelength and the wavelength range of emission light. The 
transmittance of the single long pass filter at the excitation wavelength was 2%. We used 
three long pass filters in series and obtained 0.0008% transmittance at the excitation 
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wavelength. The measured temporal profiles were the results of convolution between the 
temporal response for the ideal impulse light source and the instrument function including the 
finite pulse width of the light source and the time-resolution of the TCSPC system. Because 
the forward calculation assumes the ideal impulse source the mean time-of-flight data for the 
ideal impulse source were obtained using the following simple Eq. (15) [25,36], 



< t >=< t >, 



<t>, 



(15) 



where <t>,„ is the measured mean time-of-flight and <t>jRp is that of the instrumental function. 
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Fig. 3. Schematics of the picosecond time-resolved measurement system [35]. Long pass filter 
for measurement of emission light was removed for measurement of excitation light. 
Arrangement of the long pass filter in this figure was for mouse experiment. For phantom 
experiment the long pass filter was attached on the phantom surface. 
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Fig. 4. The transmission spectrum of the long pass filter IR-82 as well as the excitation 
wavelength and the wavelength range of emission light. 



3. 3 Phantom experiments 

Figure 5 shows a typical cylindrical phantom used in this study. The phantoms were made of 
polyacetal resin which had the scattering characteristics similar to biological tissues. The 
reduced scattering and absorption coefficients of polyacetal resin at the wavelength of 800 nm 
were fi s ' = 0.863 mm -1 and fi a = 0.0006 mm -1 respectively. The diameter and height of the 
phantoms were 30 mm and 85 mm, respectively, and the phantoms had one or two cylindrical 
target holes with a height of 65 mm and a diameter of 4 mm, 5 mm or 6 mm where 0.8% 
Intralipid solutions containing 0.25 uM to 2.0 pM ICG were injected. The 0.8% Intralipid 
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solutions had the scattering characteristics similar to polyacetal resin. Table 1 summarizes the 
conditions of 10 phantoms used in the phantom experiments, where the origin of the x-y 
coordinate was at the cylinder axis. 
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Fig. 5. Schematics of the phantoms. The long pass filter was attached on the phantom surface, 
and its position was shifted according to the measurement of the excitation or emission light. 
The long pass filters were placed on the detector plane only when measuring the emission light. 



Table 1. Summary of 10 phantoms used in the experiments 



Phantom No. 


Number of 
targets 


Diameter of target 
[mm] 


Position of target center 
(x,y) [mm] 


ICG in target 
[pM] 


1 




4.0 


(-7.5, 0) 


1.0 


2 




5.0 


(-7.5, 0) 


1.0 


3 




6.0 


(-7.5, 0) 


1.0 


4 




6.0 


(-7.5, 0) 


2.0 


5 




6.0 


(-7.5, 0) 


0.50 


6 




6.0 


(-7.5, 0) 


0.25 


7 




6.0 


(-11.5,0) 


1.0 


8 


2 


6.0 


(-6.0, 0), (6.0, 0) 


1.0, 1.0 


9 


2 


6.0 


(-5.0,0), (5.0,0) 


1.0, 1.0 


10 


2 


6.0 


(-4.0, 0), (4.0, 0) 


1.0, 1.0 



Twelve fiber bundles were used for both the source and detector and the plane of the 
detector fiber bundles was 5.5 mm below that of the source fiber bundles because of the 
geometrical limitation. When measuring the emission light the long pass filter stated above 
was put on the phantom surface as shown in Fig. 5. 

The image reconstruction was carried out with the homogeneous background reduced 
scattering coefficient, pj = 0.863 mm -1 , which was held constant during the inversion process. 
It started with the initial guess of the homogeneous absorption coefficient, [i a = 0.0006 mm -1 . 
The update was iterated until the reconstructed values xt was found to converge by checking 
that the norm of dxk was sufficiently small. 

3. 4 In vivo mouse experiment 

Figure 6 shows the schematics of an in vivo mouse experiment. In order to measure the weak 
fluorescence light it was necessary to keep good contact between the fiber bundles and the 
mouse skin. But, the mouse skin was very soft and it was difficult to keep good contact during 
the whole experimental period. In addition it was very difficult to attach a total of 16 fiber 
bundles in one plane around the deformable mouse abdomen with a diameter less than 30 mm. 
Then, an anesthetized mouse was fixed rightly in a mouse holder of an annular black 
polyacetal cylinder with inner and outer diameters of 23 mm and 30 mm, respectively, and 
with a length of 40 mm. A total of 16 fiber bundles were attached onto the surface of the 
mouse holder with an equal spacing of 22.5° in one plane perpendicular to the cylinder axis as 
shown in Fig. 7(a). In the image reconstruction process, the whole volume consisting of the 
mouse and mouse holder was the object for fluorescence tomography. To avoid light 
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propagation along the surface of the mouse holder due to bad contact between the fiber 
bundles and the mouse holder, a black soft rubber covered the whole surface of the mouse 
holder except the positions of the fiber bundles. 

The fluorophore target was a capsule containing 1 uM ICG mixed in a 0.8% Intralipid 
solution. The plexiglass capsule had inner and outer diameters of 4 mm and 6 mm, 
respectively, and a length of 30 mm as shown in Fig. 7(b). We intended to simulate a grown 
tumor accumulating a fluorophore by using the ICG capsule. The experiment was carefully 
carried out to place the capsule rightly at the position as illustrated in Figs. 6 and 7, but the 
positional error might have occurred during the preparation. We ignored the influence of the 
plexiglass because it was thin and strong diffusion effect by the tissues must overwhelm the 
influence, but there is a possibility that the light propagation was slightly influenced by the 
non-scattering and non-absorbing plexiglass resulting in slightly stronger emission light and 
higher fluorophore concentration of the reconstructed target. A long pass filter to cut the 
excitation light when measuring fluorescence light was installed in a switching device of the 
detector fiber. 



Picosecond time-resolved 
measuring sys tem 




Fig. 6. Experimental setup for mouse experiments. 

The mouse experiment was conducted with the approval of the committee of Tokyo 
Metropolitan Institute of Medical Science and in accordance with the institute guidelines for 
animal experiments. A mouse (JclTCR, 17 week old, 43.0 g) was weakly anesthetized by 
breathing diethyl ether first, and deeply anesthetized by injecting pentobarbital into the 
abdomen. The hair around the abdomen was removed by a hair removing cream. The 
fluorophore capsule was implanted beneath the skin with its center located at the position of 
(x, y, z) = (0, 8, 0) when the center of the whole object was the coordinate origin, (x, y, z) = (0, 
0, 0). The whole object including the mouse and mouse holder attached with the fiber bundles 
was fixed in a dark room to avoid stray light and time-resolved measurement was performed. 
The accumulation time for measuring one temporal profile was 15 sec, and it took about 45 
min to acquire a set of the temporal profiles necessary for reconstructing one image, because 
only one PMT was used for detection and the 16 detector fibers were switched manually. By 
some modification of the system for simultaneous use of 16 PMTs, it is possible to reduce the 
acquisition time to a few minutes which are short enough for practical use. 

As in the case of the phantom experiments, the image reconstruction was carried out with 
the homogeneous background reduced scattering coefficient, fi s ' = 0.863 mm -1 , which was 
held constant during the inversion process. It started with the initial guesses of fi a = 0.003 
mm -1 and 0.0006 mm -1 for the mouse and mouse holder, respectively. 

One updating process took about 4 seconds, and the updating process was terminated 
when the difference in the residual errors between two successive iterations was smaller than 
0.003% of the initial residual error. 
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(a) (b ) (c) (d) 



Fig. 7. (a) Picture of an anesthetized mouse fixed in a mouse holder, (b) size of the mouse 
holder and position of the ICG capsule embedded in the mouse abdomen with the measurement 
plane indicated by a vertical line, (c) arrangement of the source and detector fiber bundles; 8 
source fiber bundles (SI to S8) and 8 detector fiber bundles (Dl to D8) were alternatively 
attached onto the surface of the mouse holder with an equal spacing, (d) size of the ICG 
capsule containing 1 uM ICG mixed in a Intralipid solution. 

4. Experimental results and discussions 

4. 1 Phantom experiments 

4. 1 . 1 Measured temporal profiles 

Figure 8 shows the measured temporal profdes of the excitation and fluorescence emission, 
T x (r h ,t) and T,„(r b ,t), as well as those of the total light and zero-lifetime emission, Tj{r b ,t) and 
r*,„(r b ,0, for the case of phantom No. 3. Figures 8(a) and 8(b) are the cases of short and long 
distances between the source and detector positions, respectively. The measured excitation 
light, T x (rb,t), in Fig. 8(b) varies more slowly and its peak shifts to later time than that in Fig. 
8(a) because it is more strongly affected by scattering. A similar behavior is observed for the 
measured temporal profdes of fluorescence emission, but the values of the fluorescence 
emission are about one-order of magnitude smaller than those of the excitation. The temporal 
profdes of r*,„(rb,0 are accompanied with noises due to the deconvolution process which 
enhanced noises, but the noise levels were acceptable for image reconstruction. 



(a)detector 
sourc( 




(a) (b) 

Fig. 8. Measured temporal profiles of the excitation and emission light, r v (r b ,f) and Y m (r b ,t), 
respectively, together with the calculated zero-lifetime emission and total light, T*,„(r b ,t) and 
Tj(r b ,t), respectively, for the case of phantom No. 3. (a) and (b) are for the cases of the 
detectors close to and far from the source, respectively. 

In our previous report [21], we investigated the effect of the noise. The noise gave no 
effect on the reconstructed images for S/N ratio smaller than 40 dB. When S/N ratio was 25 
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dB, the target was reconstructed, although the absorption coefficient was reconstructed less 
than 40% of the true value. 

4. 1 .2 Results of image reconstruction 

Figure 9 shows the reconstructed images of fluorophore concentrations in [uM] for 10 
phantoms in Table 1. The color bars are normalized by the original concentration, and the 
graphs show the concentration profdes along the x axis with the red lines showing the true 
profdes. 

Figures 9(a) to 9(c) compare the effect of the target size. Although the positions of the 
target can be identified correctly the reconstructed peak concentrations increase with the 
target size, i.e., the peak concentration for phantom No. 1 with a target diameter of 4 mm is 
about 25% of the true value while that for phantom No. 3 with a target diameter of 6 mm is 
about 60% of the true value. This trend is consistent with the results of simulation in the 
previous report [21]. 
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Fig. 9. The reconstructed images of ICG concentration for the ten phantoms listed in Table 1. 
The unit of the color bars is [uM]. 

Figures 9(c), 9(d), 9(e) and 9(f) shows the results for phantoms No. 3, 4, 5 and 6, 
respectively, illustrating the effect of the target concentration with the fixed size and position 
of the targets. The position of the target was correctly identified in each case. The 
reconstructed peak concentrations are about 60% of the true values for phantoms No. 3 (true: 
1.0 uM) and No. 4 (true: 2.0 uM) and about 80% for phantoms No. 5 (true: 0.5 uM) and No. 6 
(true: 0.25 uM). These results indicate that the reconstructed peak concentration depends on 
the true concentration, but the previous simulation results did not show the dependency [21]. 
This disagreement between the phantom experiments and simulation results may be caused by 
a wrong estimation of the specific extinction coefficient, e. The concentrations are calculated 
by dividing the subtraction between fi a {r) + e N{r) and fi a {r) by e. In this study e was given 
constant e = 0.04 mnT'uM -1 for all the cases according to the literature [32] which also 
reported the dependency of e on the ICG concentration over 1 uM due to self-absorption. It 
will be necessary to correctly measure not only e, but also other fluorescence characteristics of 
lifetime, r, and quantum yield, y, which are important constants to calculate the total light <£>, 



#178955 - $15.00 USD Received 31 Oct 2012; revised 12 Jan 2013; accepted 14 Feb 2013; published 28 Mar 2013 
(C) 2013 OSA 1 April 2013 | Vol. 4, No. 4 | DOI:10.1364/J3OE.4.000635 | BIOMEDICAL OPTICS EXPRESS 647 



from Eqs. (4) and (5), toward more precise validation of the total light approach for 
fluorescence tomography. 

Figure 9(g) shows the image for phantom No. 7 in which the position of the target was 4.0 
mm closer to the surface than phantom No. 3. The reconstructed peak concentration is about 
90% of the true value, which is better than that for phantom No. 3. But, an artifact with a 
negative concentration appeared in the central area of the image. Generally speaking, 
fluctuations in the reconstructed distributions are often observed in images reconstructed by 
inversion processes, and it will be necessary to employ some regularization method to 
suppress the fluctuation and artifacts [37]. Another reason of the artifacts may be the 
employment of the photon diffusion equation which is the approximation of the radiative 
transfer equation. The photon diffusion approximation is valid for the area and time where 
photons have propagated more than a few millimeters. Therefore, when the target exists very 
close to the surface the photon diffusion equation may not be an appropriate propagation 
model [27]. 

Figures 9(h), 9(i) and 9(j) show the reconstructed images for phantoms No. 8, 9 and 10 
having two targets with different center-to-center separation (CCS) of 12 mm, 10 mm and 8 
mm, respectively, but the same size and ICG concentration. The two targets are clearly 
separated and the reconstructed peak concentration was about 60% of the true value in the 
image of phantom No. 8 shown in Fig. 9(h) while the two targets were not clearly separated in 
Figs. 9(i) and 9(j) with the CCS of 10 mm and 8 mm, respectively. In the previous simulation 
[21] two targets were clearly separated even when the CCS was as small as 6 mm although 
not shown in this paper. The degraded spatial resolution in the experiment must be due to 
larger noises in the measured temporal profiles than those assumed in the simulation. The 
spatial resolution may be improved by employing a regularization method [37] as well as the 
variable regularization parameter depending on the depth from the surface [38]. 

4.2 In vivo mouse experiment 

4.2.1 Measured temporal profiles 

The arrangement of the source and detector fiber bundles is shown in Fig. 10(a) where eight 
sources and detectors, SI to S8 and Dl to D8, respectively, were attached onto the outer 
surface in the central plane of the mouse holder alternatively with an equal spacing. The 
temporal profiles of the measured powers of the excitation and fluorescence emission light as 
well as those of the zero-lifetime emission and total light are shown in Fig. 10(b) to Fig. 1 1(e). 
The ordinates are the photon counts measured by the TCSPC system, and only Fig. 10(c) for 
fluorescence emission is in a linear scale because the photon counts were much smaller than 
other light which are in logarithmic scales. Each figure shows the cases where (i) both the 
source and detector were close to the target (S8-D6), (ii) the source was close to the target 
while the detector was far from the target (S7-D3) and vice versa (S3-D7), and (iii) both the 
source and detector were far from the target (S4-D2). When the source or detector was far 
from the target the magnitudes of the fluorescence emission were in the same order of noises. 
Resultantly, the noises of the zero-lifetime emission and total light were observed larger than 
those of the excitation light. Note that the temporal profiles of the zero-lifetime emission and 
total light were smoothed due to Wiener filtering. The integration period of time for 
calculating the mean time-of-flight of the excitation light, <t x >, in Eq. (12) was determined to 
be from the time origin to the time when the temporal curve decreased to the noise level after 
the peak. The same integration period of time was used for calculating the mean time of flight 
of the total light, <t T >. 

When looking at the temporal profiles of fluorescence emission, Fig. 10(c), the magnitude 
of the profile for the case of (ii) (S7-D3 and S3-D7, the source was close to the target while 
the detector was far from the target and vice versa) was almost the order of noises. Then, for 
the case of (iii) (S4-D2, both the source and detector were far from the target) the fluorescence 
emission from ICG would not be observed. But, the fluorescence emission of the case (iii) 
was stronger than that of the case (ii). The emission of the case (iii) must be autofluorescence 
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from the mouse tissues. This conjecture is proved by the fact that very similar strong 
emissions were observed for other combinations of the source and detector with the same 
distance as that of S4-D2, i.e., other combinations are S1-D2, S2-D3, S3-D4, S4-D5, S2-D8, 
S3-D1, S5-D3, S6-D4. There was a possibility that the light intensity observed by the 
combination of S4-D2 was a leakage of the excitation light. But that light was not eliminated 
by increasing the attenuation of the optical filter at the excitation wavelength. This fact 
indicated that the wavelength of the light was different from that of the excitation light. So the 
light observed by S4-D2 must be auto fluorescence. Image reconstruction was performed with 
the emission data including autofluorescence. 
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Fig. 10. (a) The arrangement of the source and detector fiber bundles, and the measured 
temporal profiles of (b) the excitation, (c) fluorescence emission, (d) zero lifetime emission and 
(e) total light. 

4.2.2 Reconstructed image 

The reconstructed ICG concentration distribution is shown in Fig. 1 1 where large and small 
white circles indicate the interface between the mouse and mouse holder and the true position 
of the ICG capsule, respectively. In the reconstruction process, the absorption coefficient of 
the mouse holder was fixed as that of polyacetal resin, and the ICG concentration within the 
mouse holder was fixed as zero. An area with high ICG concentration was reconstructed 
almost right on the true position of the ICG capsule with the maximum concentration of 0.94 
uM, the area occupying the concentration higher than the half maximum of about 1.6 times 
the true area, the integrated ICG concentration over this area of about 1 .2 times the true value, 
and the position shift of the maximum concentration of 2.5 mm from the true position. 




X [mm] 

Fig. 11. Reconstructed image of ICG concentration for mouse experiment. The unit of the color 
bar is [uM]. 
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Artifacts with the maximum concentration of 0.91 uM appeared along the interface 
between the mouse and mouse holder. It is known that much of the autofluorescence in 
animals is mainly emitted from the skin and colon [39]. Therefore, the major reason of the 
artifacts in Fig. 1 1 must be autofluorescence from the mouse skin because the colon in the 
mouse was out of the measuring plane. Other possible reasons of the artifacts are irregular 
contact between the fiber bundles and outer surface of the mouse holder and possible small 
gaps between the mouse skin and inner surface of the mouse holder. These irregular contacts 
will bring irregular light propagation which lead to deformation of the temporal profiles and 
to artifacts. In the mouse experiment these irregular contacts were eliminated as much as we 
could by careful setup of the experimental. Although artifacts appeared in the vicinity at the 
interface between the mouse and mouse holder, the ICG target was reconstructed reasonably 
well as in the phantom experiments. 

In the phantom experiments, we employed polyacetal resin for the background phantom 
material which has very small absorption, /u a = 0.0006 mm -1 , while the additional absorption 
by ICG was sN = 0.01 mm" 1 with e = 0.004 mm"VM"' and N = 0.25 uM. Therefore, the 
contrast, i.e., the ratio of the absorption coefficient of the target to that of the background, was 
as large as 16.7 which may be easy to be reconstructed. But in vivo experiments using small 
animals, the background absorption coefficient is much larger in the order of 0.001 mm -1 
which leads to low contrast targets. Larger absorption by the background tissues in the mouse 
experiment than in the phantom experiments seems to degrade the quality of the reconstructed 
images, because the S/N ratio of the measured data was lowered. The reconstructed target, 
however, was well localized and had the ICG concentration higher than the true value in 
contrast to the results of phantom experiment. This can be another effect of the background 
with strong absorption. When the background has a large absorption coefficient, the light 
propagation path from the source to the detector becomes narrower. Narrower path makes it 
easy to identify the position of the target. As a result, the position and area of the target can be 
localized well. And it has been shown in the literature [37] that when the reconstructed target 
is well-localized, the absorption coefficient becomes larger. Another reason of the 
reconstructed ICG concentration higher than the true value may be the existence of the 
plexiglass capsule as stated above. 

On the other hand, the small quantum yield of ICG in Intralipid (y = 0.03) may become 
larger when ICG is combined with protein in blood, for example with high density lipoprotein 
[33]. This will make the fluorescence light stronger and the image reconstruction easy in small 
animals even with highly absorbing background. New fluorophores with high quantum yields 
in the near-infrared wavelength range are actively being developed for experiments using 
small animals. Employment of new fluorophores will improve the performances of 
fluorescence tomography in general. 

5. Conclusions 

Experiments using phantoms and mouse with the embedded ICG targets were conducted to 
demonstrate the total light approach for time -domain fluorescence tomography, and the 
fallowings were concluded. 

The temporal profiles of the zero-lifetime emission and total light were obtained from the 
measured temporal profiles of the excitation and fluorescence emission light with removing 
the effect of the instrument function. The tomographic images of the ICG concentration were 
obtained from the measured mean time-of-flight data of the excitation and total light. In the 
reconstructed images of the phantom experiments, the positions, sizes and ICG concentrations 
of the embedded targets were reconstructed with the qualities similar to those in the previous 
simulation although the spatial resolution was worse than that of the previous simulation. In 
the reconstructed image of in vivo mouse experiments, the ICG target was again reconstructed 
reasonably well. However, the autofluorescence from the mouse body produced the artifacts 
in the vicinity of the mouse skin, and it is suggested that small gap between the mouse and 
inner surface of the mouse holder may lead to artifacts. 
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The total light approach for time domain FDOT was experimentally demonstrated, and 
further studies are expected for reconstruction of multiple fluorescence characteristics with 
full use of time domain data. 
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